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Abstract. By considering the scaling behaviour of various Feynman graphs 
at leading order in large Nt at the non-trivial fixed point of the (i-dimensional 
/5-function of QCD we deduce the critical exponents corresponding to the 
quark, gluon and ghost anomalous dimensions as well as the anomalous 
dimensions of the quark-quark-gluon and ghost-ghost-gluon vertices in the 
Landau gauge. As the exponents encode all orders information on the per- 
turbation series of the corresponding renormalization group functions we 
find agreement with the known three loop structure and, moreover, we pro- 
vide new information at all subsequent orders. 
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One of the more difficult exercises in understanding the quantum field 
theories, such as quantum chromo dynamics, which underpin the dynamics 
of the particles we observe in nature, is in gaining knowledge of their per- 
turbative structure to very high orders. For example, following the early 
one and two loop work of [jj], ||, [| the renormalization group functions of 
QCD are only known to third order, ||. To a large degree one is hindered 
from pushing such calculations to a further order, which will be necessary 
due to the increase in experimental energies envisaged at LHC and the SSC, 
by the huge numbers of Feynman diagrams one would have to compute, a 
subclass of which possess an intricate integral structure. Indeed such calcu- 
lations can only be tackled by ingenious algorithms which are executed on 
computer. Whilst the basic functions of the renormalization group equation 
are ordinarily calculated at successive orders in the perturbative coupling 
constant, which we will denote by g, it is not the only parameter with which 
one can compute with. For theories, such as QCD which possess Nf flavours 
of quarks one can also expand the functions in powers of 1/Nf and com- 
pute various coefficients. The technique to achieve this was introduced in 
[^, for the O(N) bosonic a model and later extended to the a model on 
CP(N), which is a U(l) gauge theory, and then to QED, || [TQ| . Basically, 
it entails computing various critical exponents at the d-dimensional fixed 
point of the theory, defined as the non-trivial zero of the /5-function, where 
the theory has a conformal symmetry, ||. The exponents which emerge 
depend on the space-time dimension, d, and are computed order by order 
in powers of 1/-/V/. Further, as they are calculated at a fixed point the 
renormalization group equation takes a simpler form there to the extent 
that exponents have simple relations with the corresponding renormaliza- 
tion group equations, such as the wave function or vertex renormalization, 
at criticality, [pTlj] . Moreover, one can undo such relations to deduce infor- 
mation on the perturbative structure at non-critical g. In particular the 
coefficients of each power of the coupling constant, at the order in large N 
one is working to, can be determined and they correspond to those of the 
MS scheme. This can be understood from the fact that at the critical point 
there is a conformal symmetry and thus the exponents must correspond to 
a mass independent scheme. Another success of the method of [||, fjj] is the 
possibility of pushing the calculations to 0(1/Nj), flotl , which is a depth 
not possible by conventional large N methods. 

The purpose of this letter is to set up the analogous approach for QCD, 
extending the earlier QED work of M. The aim is to deduce non-trivial 
exponents for various fields in QCD and demonstrate their agreement with 
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the known three loop perturbative MS calculations, [Q, [|, and then to deduce 
information at higher orders. By non-trivial we mean determining exponents 
which are not merely deduced by examining the analogous QED results and 
multiplying those by the appropriate non-abelian group factors. Such non- 
trivial results will be derived either from graphs involving ghost fields, the 
triple gluon vertex or those which would otherwise be absent in QED by 
Furry's theorem. In particular, we will deduce the anomalous dimensions 
of the quark, gluon and ghost fields at leading order in large Nf as well 
as the vertex anomalous dimensions of the quark-quark-gluon (qqg) and 
ghost-ghost-gluon (gcc) vertices, in the Landau gauge. Although these field 
anomalous dimensions are gauge dependent they must be computed initially 
before attempting to determine gauge independent quantities such as the (5- 
function which will be a considerably more involved computation and is our 
long term goal. This is akin to the perturbative approach when the wave 
function renormalization of each field must be determined first. 

The version of the QCD lagrangian we use to compute exponents is 

(F a ) 2 - 1 

L = — - J ^r~ + i^W 1 + A^Tfyf^ - -^f ahc d ll A a v A^ b A vc 

- d^c a d^c a + f abc d^c a c b A» c (1) 

where if) 11 is the quark field, l<i<Nf,l<I< N c , A a ^ is the gluon field, 
1 < a < AT? - 1, F^, = d^Al - d v A°, c a and c° are the ghost fields, f 
is the covariant gauge parameter, Tfj are the generators of the non-abelian 
gauge group, f abc its structure constants and e is the coupling constant. It 
is important to note that we have absorbed a power of the coupling constant 
into the definition of the gauge field so that the qqg vertex has the same 
style of interaction as the electron photon vertex of QED, 0, [n|, which was 
primarily for applying the techniques of uniqueness of [[j^] to computing 
massless Feynman diagrams which arise at 0(1/Nj). To make more explicit 
the notions already discussed we note that the fixed point where we will 
concentrate our analysis is given by [3{g c ) = 0, where g c is the non-trivial 
zero of the d = 4 — 2e dimensional /3-function. It has been calculated in MS 
using dimensional regularization and in terms of the dimensionless coupling 
constant g = {e/2ir) 2 , in the notation of [13], it is, [1-5], 



(3(g) = {d-4)g + 



2 -T{R)N f - ^C 2 (G) 



g 2 
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+ 



l -C 2 (R)T{R)N f + ^C 2 (G)T(R)N f - ^C 2 2 (G) 
^-C 2 (R)T\R)Nf + ^-C 2 (G)T\R)Nf 



72 

+ ±C 2 2 (R)T(R)N f 



432 
205 



C 2 {R)C 2 {G)T{R)N f 



1415 



2857 



864 

from which it follows that 
3e 1 



C$(G)T(R)N f + ^CKG) 



g 4 + o(g 5 ) 



9c 



+ 



T(R)N f AT 2 (R)N 2 



33C 2 (G)e - (27C 2 (i?) + 45C 2 (G)) e 2 



+ 



+ o 



N f, 



(3) 



The Casimirs for a general classical Lie group which appear in (2) and (3) 
are defined via 



Tr(T a T b ) = T(R)5 ab , T a T a = C 2 (R)I , f acd f 



tab 



acd rbcd 



C 2 {G)5 ah (4) 



and for SU(N C ), T(R) = |, C 2 {R) = {N 2 - l)/2N c and C 2 {G) = N c . 

At g c the fields of (1) obey simple power law behaviour in analogy with 
ideas in statistical mechanics and therefore in momentum space the asymp- 
totic scaling forms of the propagators of each field in the critical region as 
k 2 — > oo are, ||, 



c(k) 



Aft 



(k 2 y- a 
c 



B 



(k 2 ) 



Via* 



kfik u 



k 2 



(k 2 ) 



M-7 



(5) 



The quantities A, B and C are the respective fc-independent amplitudes 
of each field and a, (5 and 7 are their dimensions which are built out of 
the canonical (or classical) dimension and the anomalous contribution. The 
latter arises as a result of quantum corrections such as radiative corrections 
to the underlying Green's functions. In particular, we define the exponents 
to be 



a 



(J,-l + kv , = l-i]-x , 7 



(JL—l + ^Tjc 



(6) 
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where d = 2/i and 77, r\ c and % are respectively the quark, ghost and qqg 
vertex anomalous dimensions. Each is 0(1/Nf) and depends on /i and N c 
through the Casimirs of (4) and we will determine their leading order large 
Nf structure here. Also their relation to the conventional renormalization 



group functions are, [13], 



V < — > 72 (s c ) , V + X < — ► l3(g c ) 

X < > 71 (9c) , Vc < > 73 (&) (7) 

where 72(5) is the quark wave function renormalization, 71(5) is the quark 
gluon vertex renormalization function, 73(5) is the gluon wave function 
renormalization and 73(5) corresponds to the ghost wave function renor- 
malization. The third relation of (6) arises from analysing the ghost kinetic 
term in (1) but it can be related to other anomalous dimensions by exam- 
ining the ghost interaction term to obtain 

Vc = V + X ~ Xc (8) 

where Xc is the gcc vertex anomalous dimension being defined in a similar 
way to x- This result, which we will use later to deduce r\ c and demonstrate 
its agreement with perturbation theory, is nothing more than an expression 



of one of the Slavnov Taylor identities for QCD in exponent language, [14]. 
A similar observation concerning the QED Ward identity was made in [15]. 
As we are working in a gauge theory we must choose a particular gauge to 
calculate within. As the covariant gauge parameter remains unrenormalized 
in the Landau gauge, we have made the choice £ = 0, ||] . This is important 
for comparing the structure of gauge dependent exponents with perturbation 
theory which we do later. 

The method to solve for (7) involves examining the Dyson equations of 
various Green's functions truncated to 0(1/Nt) in the critical region where 
the propagators of the Feynman graphs are replaced by the expressions (5). 
Although the analysis of the quark and gluon 2-point functions is trivial in 
the sense defined earlier, it is worth discussing it to illustrate the simplicity 
of the method and the Dyson equations are shown in fig. 1 where dotted lines 
correspond to quarks and solid lines to ghosts. We use dressed propagators 
since the effect of including a non-zero anomalous dimension in (5) is to 
reproduce the infinite chain of bubble graphs one would ordinarily have to 
consider, which was the original approach to examining QED in the large 
Nf expansion, Using (5) means fewer Feynman graphs have to be 

analysed. Several points concerning the graphs of fig. 1 ought to be noted. 



5 



First, there are no graphs in the gluon 2-point function involving ghost loops, 
triple or quartic gluon vertices since these are 0(1/Nf) with respect to the 
quark loop. (A similar set of graphs were analysed in a different context 



in [ IT | . ) Second, the quantities tp , A~^ and c _1 denote the respective 
2-point functions and their asymptotic scaling forms can be deduced from 
(5) by simple momentum space inversion where the gluon is inverted on the 
transverse subspace, M, as 



kfxk v 



A(k 2 ) a ~t J,+1 ' ^ v ' B{k 2 )P~ 

c~\k) ~ — =r (9) 

Thus the first two graphs of fig. 1 can be represented in the critical region 
by 

= 1 + ~ l M± l ~J* )a a{^ - t3)a(ji - a)a(a + /3)C 2 (R)z (10) 

8a 3 T(R)Nf „. . , 
= 1+ /n / fl 2 /t-aa2a + lz 11 
(2a + 1) 

where only the transverse part of the gluon equation is physically relevant, 
0, trl = 4, z = A 2 B (and y = BC 2 ) and we have set a(a) = T(fi — a)/T(a) 
for all a. In (10) and (11) the only unknowns at leading order in 1/Nt are 
rji and zi, where for example rj = Ya^Li Vi/Nf with rji = rji(fx, N c ). Although 
it may appear that xi a l so arises through the factor a(a + (3) in (10), in the 
momentum space approach to solving the Dyson equations in this critical 
point method, one sets xi = at leading order so that /3 = 1 — rj here. 
This is partly to ensure consistency between the momentum and coordinate 
space solutions for 771 , which can easily be verified by examining the x-space 
QED solution, [ficfl , or considering the same calculation in the toy O(N) a 
model of JOI]. More importantly, though, since the method we are using 
relies very much on the conformal nature of the fixed point, the qqg 3-point 
vertex is in fact only conformal or unique, [10|, at leading order when xi 



0, H 0, |To|, which justifies this point. Further, the factors (p 2 )~ x have been 
omitted from the second terms of (10) and (11) partly for this reason but 
also because xi = 0(1/Nf) and so its 1/Nf expansion, ie 1 — (xi lnp 2 )/Nf, 
means it will only be important for deducing 772. The treatment of this can 
be seen, for example, in the QED work of [O. Thus eliminating Z\ between 
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(10) and (11) yields 



m 



(2 f x-l)(n-2)T(2 t i)C 2 (R) 



4r2(^)r( A , + i)r(2- M )T( J R) (12) 

which corresponds to the QED result of Q aside from the non-abelian factor 
of C2(R)/T(R) = 8/3 for QCD and agrees with the 3-loop perturbative 
expansion of 72 (g c ). 

To deduce the vertex anomalous dimensions one carries out a similar 
analysis using instead the truncated leading order large Nf graphs contribut- 
ing to the qqg and gcc 3-point functions which are illustrated in figs 2 and 
3 respectively. However, naively computing each individual graph with the 
leading order exponents one would discover that they are infinite and there- 
fore need to be regularized. This is achieved by shifting the gluon exponent 
by an infinitesimal amount A, /?—>/? — A, with A playing a similar role 
to the e in dimensional regularization, [10]. The aim now is, following the 
ideas of p8|, to compute each graph with the shifted exponent to observe 



+ Q + Rlnp 2 + 0(A) 



(13) 



that they have the following structure at 0(1/Nf) 

P 
A 

where P, Q and R will depend on \x and N c and p is the external momentum 
flowing through the external quark or ghost fields. The simple pole of (13) 
is absorbed into the usual vertex counterterm, [18]. This leaves a A-finite 
Green's function but the scaling behaviour would be spoiled by the presence 
of the Inp 2 term. However, general arguments, [Tq |, imply that the overall 
scaling behaviour of the qqg vertex, for instance, has to be (p 2 ) x ^ 2 and 
therefore the Inp 2 terms at each order in 1/iVf must re-sum to this form 
which allows one to determine an expression for \- 

Thus substituting (5) into the graphs of fig. 2 the first graph involves 
multiplying the QED result of [15] by the appropriate group factor, so that 
it contributes 

- C 2 (R) - ^ 



V? 
T(R) 



(14) 



to Xl where r/i = C2{R)f}\ /T(R). Useful in computing the contribution 
from the two loop graphs of fig. 2 is the result 

kf,k u (a + - n - 2)a(a - l)a(/3 - l)a(2n - a - j3 + 2) 



k {k^{{k-p) 2 Y 



£_ 

2(a- l)(/3- l)(p2) Q +^ 
2pti,Pv {n~a + l)(a + /? 



P 



(/*-/?) 



(15) 
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valid for all a and j3 and they and the second graph of fig. 2 are proportional 
to C2(G) and contribute 



to Xl- Thus, we have 



Xi 



(/i-l)C 2 (GK 
2( M - 2)T(R) 

C 2 (G) 



C 2 (R) + 



2(/z - 2) 



T(i?) 



from which we deduce the gluon anomalous dimension as 



+ Xl 



2( M -2)T(i?) 



(16) 



(17) 



(18) 



To check that (17) and (18) are in agreement with perturbation theory we 
recall that the three loop MS result of || in arbitrary covariant gauge, in 
our notation, is 



7i(0) 



£C 2 (R) + Q±$-C 2 (G) 



+ 



l lc 2 {R)C 2 {G) - f Cf (A) + ^C 2 2 (G) - 2C 2 (R)T(R)N f 



24 



C 2 (G)T(R)N f - (1 - i)C 2 (G) [ -C 2 (R) + -C 2 (G) 



+ 



6 

(1-0 2 C 2 (G) 



15 



C 2 (R) + -C 2 (G) 



+ 



-Ci(R) + Ci(R)C 2 (G) 12C(3) - — 



143 



+ C 2 {R)Cl(G) 



10559 15 



144 



^rC(3) + 7 C(3) + 



10703 \ 



864 



) cm 



- 3C%(R)T(R)N f + 2C 2 (R)C 2 (G)T(R)N f ( 6C(3) - — 



(%{G)T{R)N f f 9C(3) + ^) + fc 2 (R)T 2 (R) 



Ni 



35 
27 



C 2 {G)T\R)Nf - (1 - C 2 {R)Cj(G) -C(3) + 



371 \ 
"32"; 



+ (18C(3) + 127) - C 2 (G)T(R)Nf (^C 2 (R) + C 2 (G)) 
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+ (1 - A Ci(G) (C 2 (R) £C(3) + ™j 



69 



C 2 (G) 



- (1 



32 



32 

3 



9_ 

32 



which with (3) gives, in the Landau gauge, 

C 2 (G) 



7i (5c 



T(R)Nf 

C 2 (R) 
T(R)Nf 



9 

-e 



9 



15 n 35 o 4 
16 £ "32 £ +0(£ ) 



15 



+ o 



(3C(3) + 27) 



(19) 




(20) 



Expanding (17) in powers of e = 2 — fi to 0(e 3 ), it is easy to see that it is in 
complete agreement with (20) which is a stringent check on our results. More 
importantly, the higher 0(1/NA coefficients of (19) can now be deduced by 
extending this argument. If we write the leading order part of (19) as 



7i(5) 



\C 2 {G)g + Y,an(T{R)N f Y- 1 g n 



(21) 



then (17) implies 



04 = — 



1 

36 



35 



36 



n=2 



83 \ 



-C 2 (R) + C(3)- — )C 2 (G) 



144/ 



a 5 



1 
108 



C 2 (R) (2C(3) - I 



C 2 (G) 



3C(4) 



5C(3) 



65 
48 



[22) 



It is interesting to note that the first appearance of the £ parameter in (19) 
and also in the other three loop renormalization group functions of || is at 
0{1/Nj). In fact this is true to all orders so that the leading order large 
Nf coefficients are gauge independent. Thus the exponents we have written 
down encode gauge independent information, at all orders, of part of the 
renormalization group functions. 

Repeating the analysis for the gcc vertex graphs of fig. 3, one finds that 
each is A-finite and therefore at 0(1/ N A 



Xcl 







which implies from (8) that 



Vcl 



C2(G) V f 
2Gu - 2)T(R) 



(23) 



(24) 
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Expanding (24) in powers of e as before and comparing with the two loop 
arbitrary gauge calculation of || one again finds agreement. The three 
loop ghost anomalous dimension has been calculated in Although that 
calculation was carried out in the Feynman gauge the three loop 0(l/Nf) 
coefficient is in agreement with the 0(e 3 ) term of the expansion of (24) in 
the Landau gauge and our observation concerning the non-appearance of £ 
at third order at 0(1/Nf) in previous renormalization group functions also 
hold at this order. 

Finally, it is worth recording the leading order form of a gauge inde- 
pendent exponent which encodes all orders coefficients on the perturbative 
structure of the renormalization of the quark mass, Jm(g)- It is trivial in 
the sense we mentioned earlier and is deduced by computing the anoma- 
lous dimension of the composite operator ipip within the Green's function 
(ipl^ij;]^) . This was discussed originally in for QED and including the 
structure which arises from the non-abelian nature of (1), we find that the 
exponent is 

2C 2 (R)r,° 

lM = ~ P)rp (25) 

and its e-expansion is in agreement with the MS perturbation series of [19]. 

We conclude by remarking that we have provided the first stage in the 
analysis of QCD in the large Nf expansion where the aim is to deduce the 
coefficients of the perturbative renormalization group functions in order to 
provide new results as well as checking existing results. It ought now to 
be possible to deduce the /3-function in the same approximation as well as 
building on the 0(\/Nf) QED work of @ to gain a deeper insight into the 
theory. 
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Figure Captions. 

Fig. 1. Leading order skeleton Dyson equations where dotted lines are 
quarks and solid lines ghosts. 

Fig. 2. Leading order graphs for A^^T""^ . 
Fig. 3. Leading order graphs for f abc A a d^c b c c . 
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